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Abstract 

In the present paper we consider a theory of gravity in which not only curvature 
but also torsion is explicitly present in the Lagrangian, both with their own coupling 
constant. In particular, we discuss the couplings to Dirac fields and spin fluids: in 
the case of Dirac fields, we discuss how in our approach, the Dirac self-interactions 
depend on the coupling constant as a parameter that may even make these non- 
linearities manifest at subatomic scales, showing different applications according to 
the value of the parameter we have assigned; in the case of spin fluids, we discuss 
FLRW cosmological models arising from the proposed theory. 
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^ 1 Introduction 

The issue about the torsional completion of gravity has already covered several decades, 
and still it is not fully resolved: although torsion is a natural part of the most complete 
metric-compatible connection, nevertheless there are also reasons against torsion that are 
jj] ' admittedly justified. In this paper, we are going to bring new insights in this topic. 

As we have just mentioned, once the covariant derivative is written in the most general 
way, it is defined upon a connection that in general is not symmetric in the two lower 
indices, and therefore torsion does not generally vanish; although the implementation 
of the principle of equivalence and causality may seem to impose the symmetry of the 
connection, and thus restricting torsion to be equal to zero [TJ [5] , nevertheless a deeper 
analysis has shown that even in the case in which these principles hold then torsion is 
only restricted to be completely antisymmetric [3J HI [5J |H1 the inclusion of torsion 
beside the curvature makes it possible to introduce a torsion-spin coupling beside the 
curvature-energy coupling, fully realizing the geometry-matter coupling prescription. It is 
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actually quite intriguing that if torsion is completely antisymmetric this complete coupling 
prescription could only be realized by the one fermion field that we have ever been able 
to observe at the moment 8, 9 ; and if we now focus on this single fermion field, namely 
the Dirac spinor, the torsional effects manifest themselves as self-interactions, capable of 
giving a dynamical explanation of the exclusion principle |10| . 

However, there are also problems that seem to push against the acceptation of torsion: 
the most important of which is the fact that the torsional effects influence the dynamics 
only at the Planck scale and not much beyond, making torsion negligible in almost every 
situation; this is due to the fact that torsion has a coupling constant that is the gravita- 
tional constant. The fact that torsion and curvature share the coupling constant might 
have been thought to have something to do with the fact that both fields are geometrical 
fields and therefore with the same strength; on the other hand however, despite torsion 
and curvature come in fact from a common geometrical background, nevertheless they 
are independent fields, and consequently endowing them with an identical strength is not 
justifiable a priori, and they should in general have different coupling constants. In the 
present paper, we shall propound what we believe to be one of the simplest theories of 
gravity where torsion is present with its own coupling constant, showing that the presence 
of torsion with a properly tuned coupling gives consequences already at subatomic scales. 

This theory of gravity is simply based on the idea that, as for a complete description 
of the underlying geometrical properties not only curvature but also torsion is to be con- 
sidered, then in order to study the dynamical properties the Lagrangian should not only 
be written in terms of the curvature but also torsional terms: this extension is natural and 
in literature it has already been put forward [TTJ Q21 [T21 [TH [TS]; in those papers, when all 
results are accounted, one may see that there are several terms that need be added to the 
Lagrangian, but if one assumes the least-order derivative hypothesis, then there can only 
be terms that in the curvature are linear and in torsion are quadratic, so that if one addi- 
tionally wants to neglect parity- violating terms then one is left with the usual torsionless 
Lagrangian plus three torsion-squared contributions |16j . Here however we would like to 
consider only one of these contributions, since we are not interested in developing the most 
mathematically comprehensive of the models but its physical consequences; thus instead 
of dwelling in all mathematical details of the most general theory, we shall focus on a 
specific model so to have a clearer look at the modifications that such a generalization will 
bring into the phenomenological aspects. As we will see, the presence of this additional 
quadratic-torsion term will have the effect of producing a shift in the value of the constant 
with which torsion couples to matter fields, so that the torsional coupling constant is not 
necessarily the Newton gravitational constant, and its value still to be determined may 
be large indeed, and in particular one may employ this modification to see that torsional 
effects may be present in particle physics and amplified in cosmology. 

In particular, for particle physics like in [TB] we are here going to study the seminal case 
of Dirac fields, while for the cosmological applications differently from [HI] where ELKO 
have been consider here we are going to employ the spin fluid; in addition we will employ 
a different formalism [T7J [22 [HZ] to write the structure of the theory. 

The layout of the paper is as follows: in section [2] we recall the main features of the 
^-bundles geometry, in which we develop our physical theory; in section [3] we formulate 
the square-torsion modification of the original Einstein-Cartan theory, showing that in 
this way curvature and torsion appear in the dynamical field equations each with its own 
coupling constant; in section!?] we apply this model to two of the most relevant matter 
fields with spin, the Dirac spinor and the spin fluid: for the Dirac field we discuss how, by 
choosing a suitable value for the coupling constant of torsion, it is possible to show that, 
in the case of two leptons, the Dirac equation has self-interactions that can be written 
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in the form of weak forces, while in the case of two neutrinos, those self-interactions can 
be used to model neutrino oscillation; in the case of the spin fluid, we study cosmological 
applications, showing how the already known cosmological effects of torsion (absence of 
initial singularity, initial accelerated expansion, horizon problem) may be affected by the 
choice of the torsional coupling constant. 

2 The Geometrical Framework 

To begin with, we shall recall the notation and convention we will employ in the following 
of the paper: readers that are familiar with the formalism of the ^-bundles may skip 
this sections and pass directly to the following one; those who are not, may also refer to 

03 QB Ql 1201 HQ for further details. 

Let us consider a 4-dimensional manifold with a metric tensor g of signature r\ = 
(1, 3) = (1, —1, —1, —1). Let £ be the co-frame bundle of M and P —> M be a principal 
fiber bundle over M. We denote by C := Ji(P)/SO(l, 3) the space of principal connections 
over P, referring £ and C to local coordinates x % , ef (£, fi = 1, . . . , 4) and x l , uj i ^ (fj, < v). 

The configuration space of our theory is the fibered product £xmC (£ x C for short) over 
M. The dynamical fields are (local) sections of £ x C, namely pairs formed by a tetrad 
field e(x) = ef(x)dx l and a spin connection 1-form oj(x) = uj^" (x) dx l , automatically 
metric-compatible with the metric g(x) = T]^ e^(x) ® e u (x). 

We consider the first J"-bundle J(£ x C) associated with the fibration £ x C — > M; 
it is a fiber bundle built similarly to an ordinary jet-bundle, but the first order contact 
between sections is considered with respect to exterior differential. 

The starting point is be the first jet-bundle J\(£ x C) associated with the fibration 
£ x C — ► M, referred to local jet-coordinates x t ,ef,uj i ,l ' / ,ef-,u) i - iiv . Then, we introduce 
on Ji(£ x C) the following equivalence relation: let z = {x l , ef , ui i ■, u^^") and z = 
(x\ if , £) i ^ u , efj : , Co i j iV ) be two elements of J\(£ x C), having the same projection x over 
M; denoting by (e M (x), u^ v (x)) and (e M (x), Lb^ v (x)) two different sections of the bundle 
£ x C — > M, respectively chosen among the representatives of the equivalence classes z 
and z, we say that z is equivalent to z if and only if 

e"(a;) = g"(as), u*"(x) = u>^{x) (2.1a) 

and 

de^(x) = di^(x), Duj^(x) = Dlu^(x) (2.1b) 

where D denotes the covariant differential induced by the connection. In local coordinates, 
it is easily seen that z ~ z if and only if the following identities hold 

K - 4) = K - K- M1/ - V) = ( V - ^-D- ( 2 - 2b ) 

The quotient space J±(£ xC)/ ~ is denoted by ,/(£ x C) and the corresponding canonical 
projection is indicated by p : Ji(£ X C) — s> ^7(f x C). A system of local fibered coor- 
dinates on the bundle J{£ x C) is provided by the set of functions x l ,ef,LJ i ^ v ,Efj := 

\ (efj - e^) ,^ ii A ' l/ := \ (w,/" - uj^") (i < j). The quotient projection p endows the 
bundle J{£xC) with most of the standard features of jet-bundles geometry (^-extension of 
sections, contact forms, J'-prolongation of morphisms and vector fields), which are needed 
to implement variational calculus on J(£ x C). Referring the reader to jTTJ [TSJ HH [2D] for a 
detailed discussion on J"-bundles geometry, the relevant fact we need to recall here is that 
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the components of the torsion and curvature tensors can be chosen as fiber J^-coordinates 
on J(f x C). In fact, the definitions of torsion and curvature components 

^= 2 ^ + ^Vi-^Vi (2-3a) 

R {j ^ = 2n 3 r + u> t >> - w/ x w 4 Xv (2.3b) 

can be regarded as fiber coordinate transformations on J{£ x C), allowing to refer the 
bundle J{£ x C) to local coordinates x l , e , u i ^ v , (i < j), R^ ^ v (i < j, \i < v). In 
such coordinates, local sections 7 : M — > J{£ x C) are expressed as 

7 : x -> (x* , ef (x) , Wj (x) , T£ (x) , iZy (x) ) (2.4) 

In particular, a section 7 is said holonomic if it is the j7-extension 7 = J a of a section 
<7 : M — ► £ X C. In local coordinates, a section is holonomic if it satisfies the relations 17 

de^tx) f) P ^(T\ 



r, i>n, , ( X ) ® x i ''"fa) . ,, , \,- , , , \, , , ,,. rl , 

« ^ = — d? d& — ' x\ x ) u i W ~ (x) (2.5b) 



namely if the quantities T^-(x) and i?^ ^(x) are precisely the components of the torsion 
and curvature tensors associated with the tetrad e^(x) and the connection ^ (x), in 
turn representing the section a. 

We also recall that the bundle J(£ x C) is endowed with a suitable contact bundle, 
locally spanned by the following 2-forms 

0" = de$ A dx l + dx l A dx 3 (2.6a) 

0"" = d^ A dx l + Q^" dx l A dx j (2.6b) 

It is easily seen that a section 7 : M — > J{£ x C) is holonomic if and only if it satisfies 
the condition 7*(0 M ) = 7*(6> A " y ) = V/x, u = 1, . . . ,4. Moreover, in the local coordinates 
x,e,w,T,R the 2-forms (12.61) can be expressed as 

0P = t m _ a nd 0^ = - itf 1 " (2.7) 

being r' 1 = def A dx l + u^^da? A dx 4 , T M = dx l A dx J , = duf A dx J + 
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dxJ A dx l and i?^ = ±i? ., ^ dx i A dx^'. 



A Lagrangian on x C) is any horizontal 4-form, locally expressed as 

L = £(x i ) ef ) a;.^,^,ii ii nd S (2.8) 

We can associate with any given Lagrangian a corresponding Poincare-Cartan 4-form, 
having local expression (see [17] ) 

1 BC 1 d£ 

Q = £ds- -—- 9 a A ds hk - 5 a/3 A ds** (2.9) 

where ds hk := J J ds. 

Any Lagrangian (|2.8p gives rise to an associated a variational principle of the kind 

A{a) = J Jo*(&) = J J<j*(Cds) (2.10) 
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where a : M — > £ xC denotes an arbitrary section and Jo : M — > J{£ x C) its J'-extension 
satisfying eqs. (|2.5[) . Referring again the reader to [T7] for details and comments, we just 
recall that the corresponding Euler-Lagrange equations can be expressed as 

Jo* (J(X)J dO) = (2.11) 

for all J'-prolongable vector fields X on £ x C. Due to the arbitrariness of X and the 
holonomy of the j7"-extension Jo, from the the requirement (|2.1ip we get two sets of final 
field equations 



and 



( dC dC a dC a dC a dC a 



(2.12b) 



Jo* 



3 The Gravitational Square- Torsion Lagrangian 

The geometrical setting outlined above does not only have the advantage of allowing a 
geometric meaning to general gauge theories, but also that of highlighting that, as both 
tetrad and spin connection are independent and they must be treated equally, similarly 
both torsion and curvature tensors have to be considered on an equal level of importance; 
this is certainly true during the variation of the Lagrangian that leads to the field equations, 
and it should therefore not surprise that such a complementarity has to be found in the 
construction of the Lagrangian itself. This is however not the case in general: for instance, 
the Einstein-Cartan theory, that is the simplest that can be derived within this formalism, 
comes from a Lagrangian of the form Cq — \f\g\R which clearly favours curvature tensor 
over torsion, despite the fact that both curvature and torsion are taken as dynamical 
variables in deriving the field equations. So the Einstein-Cartan theory seems too simple. 

Inasmuch as both curvature and torsion are taken as independent degrees of freedom, 
a more general treatment allowing torsion not only implicitly within the curvature but 
also explicitly on its own appears to be more elegant: an example of theories in which 
torsion must be considered explicitly and not only implicitly through the curvature so 
that the Lagrangian is constructed on an intertwined mix between torsion and curvature 
tensors are the conformal theories of gravity |22j . However, these theories have higher- 
order derivative Lagrangians, and consequently their torsion can propagate in regions in 
which the effects of torsion are expected to be vanishingly small [53]; even in the case in 
which torsion does vanish so that the field equations of Weyl gravity are recovered [231 > 
there are nevertheless solutions that do not reduce to those of the known Einstein limit 
[25l [26] and so it is not surprising that discrepancies with observations arise [27] . This is 
not a peculiarity of conformal gravity but of all higher-order theories of gravitation, and 
a wise choice might consequently be to look for theories that have least-order derivative 
Lagrangians so to have chances to recover the proper limit. 

The resulting situation is thus a balance between the issues of mathematical coherence 
and experimental consistency, the first pushing toward a generalization in which torsion 
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is included beside curvature while the second pulling to remain within the observational 
constraints, a tension resolved by a gravitational Lagrangian of the least-order derivative 
in curvature and torsion. A gravitational Lagrangian of this kind, linear in curvature and 
quadratic in torsion is 

C G = eR + ekT (3.1) 

where e = det(ef), R = R^^e^el and T = T pq a TPl. In the Lagrangian (gU), the 
fundamental constant for R is the Newton constant normalized to the unity while the 
fundamental constant for T is k yet to be determined. This gravitational Lagrangian is 
supplemented by the matter Lagrangian Cm in the full Lagrangian Cq — Cm as usual. 
Taking the identities J^r = ee^ and ^Jr = — e^e^ into account, we have 

dCr 



ee\R - 2eR*°e\ + ek {T°T™e\ + 4T*_ °T\ a e\) (3.2a) 



dC c 



2e [e*et - e^e*] (3.2b) 



- " '//' 



= Aekr,^T st °g> s g lt (3.2c) 

In view of this and supposing Cm independent of torsion and curvature, eqs. (|2.12[) become 

ee^R - 2eR^e\ + ek (T pq °T™e^ + 4T* p ^ ka e k ) + 

d , i ' \ t 8Cm ( 3 - 3a ) 



and 



d 

dx k 



dx k v ' 1 k " a de? 



> (e*et - e^e*)] + Wfc A M [2e (e^t - e* A e*)] 



fc \ [2e (e£e A - e^)] - 4 e fcT,V + 4efcT/„ = ^ 



(3.3b) 



After some calculations, eqs. p.3b[) may be rewritten in the form 



e (T ts a - T ta °e a s + TJe a t ) - 2ekT s a t + 2ekT t a s = -^-^e^ef (3.4) 



ldC_M 

2 8^ 

From p.4l) , by saturating with e l a , we have 

T ts 1 - T t 5\ + T s 5\ - 2kT s 1 1 + 2kT t \ = (3.5) 

where S^* := ^^T^e^e^ and T t := T (;5 S . By contracting the indices i and s, we obtain 

2(fc - l)T t = 5 t (3.6) 

From eq. (|3.6[) we see that the case fc = 1 is compatible only with St = 0. Supposing then 
k ^ 1, we can rewrite eqs. p. 61) in the form 

Tj - 2kT s * 1 + 2kT t \ = +—L-S t 6i - i^-^SsSl + Sj (3.7) 
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By adding and subtracting the term 2kT st i in the left-hand side of eq. Q3.7p . we get 
(1 + 2k)T tsl + 2kT lst + 2kT stl + 2kT tls = + S t g is - S s g it + S tsl (3.8) 

On the other hand, by contracting eq. (|3.7p with e % €huvz we get 

2 

2Tist + 2T st i + 2T t i S = 7- -tt (Si s t + S s ti + Su s ) (3-9) 

(1 — 4k) 

By inserting the content of (|3.9[) into eq. (|3.8p we end up with the final relation 



2k 1 

Ttsl = ~(l-4fc)(l + 2fe) (5lst + Sstl + Stls) + 2(k-l)(l + 2k) St9ls + 

-SsSit + jz — — ^rprStsi 



(3.10) 



2(fc-l)(l + 2fc) (1 + 2/c) 



where k ^ Eq. (l3~T0l) expresses the torsion tensor in terms of the spin density. 

It is worth noticing that in the case Stsi = we have zero torsion and then the theory 
reduces to Einstein's General Relativity. We stress that in the previous discussion we 
have supposed k ^ {1, j, — |}; such values of parameter k are pathological for the present 
theory and they should be avoided. Indeed, it is easily seen that for k — {1, 3, — 5}, eqs. 
(I3.5[) are unable to uniquely determine the whole torsion tensor. In detail, for k = 1 eqs. 
(I3.5P become 

Tj - T t 6l + T s Sl -2T S \ + 2T t \ = S ts * (3.11) 

and it is a straightforward matter to verify that solutions of p. lip are determined at least 
up to vector components of the form ViSj — Vj<5' 1 , for any Vj,. Analogously, for k = ^ the 
torsional equations Q3.5p are 

Tj - T t 6l + T s 8\ -\t s \ + \T t \ = Sj (3.12) 

This time, solutions of (|3.12p are determined at least up to totally antisymmetric compo- 
nents of the form Vijh = Vujh]- Finally, for k = — ^ eqs. Q3.5p assume the form 

T ts < - T t Sl + T s 5\ +T s \-T t \ = Sj (3.13) 

and admit solutions defined at least up to traceless non totally antisymmetric components 
of the kind V ljh - \ {Yi9jh - Vjg lh ) - V[ ijh j, where V l]h = -V jlh is a generic 3-tensor 
antisymmetric in the first two indices and Vi — 3 is its trace vector. 
Moreover, by saturating eqs. (|3.3ap with we get 

Rij - ~ \ {T mh T^ h 9lJ + 4T ivq T*. 9 ) + 2kV h T% + 2kK h fT% = S« (3.14) 

where, using the relationships R h ki j = ^ ' fjavef^&k anc ^ -^ij h ~ ^ij amon g the 
curvature and torsion tensors related respectively to the spin connection u>, M1/ and the 

associated linear connection T i:j h = e 1 ^ f^j- + u> i , := R h ih j and R = R % i are the 

Ricci tensor and scalar curvature associated with the linear connection T, denotes the 
Levi-Civita covariant derivative associated with the metric tensor £ l „ := — j- d Jr%' 

3 L J> J 2e ae'r J 

indicates the matter energy-impulse tensor and 

1 



Kii n = 7 i {-T ii h + T j \-T\ J ) (3.15) 
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represents the usual contorsion tensor. Recalling the decomposition of the metric compat- 
ible connection T in terms of the Levi-Civita connection T and contorsion tensor K 

r iJ h = r iJ h -K ij h (3.16) 

we can decompose the Ricci tensor as 

Rl] = R l3 + V jKhi h - V h K 3l h + K/K hp h - KjK Jp h (3.17) 

and then, inserting the content of eqs. (I3.15P and (|3.16p into eqs. (|3. 14[) . we can separate 
the Levi-Civita contributions from the torsional ones. 

To conclude this section, we recall the conservation laws holding for energy and spin 

V^ + ^-E^-i^. = (3.18) 

y h S tsh + T h S tsh + S* s - S s * = (3.19) 

The latter result from the invariance of the Lagrangian under diffeomorphisms and Lorentz 
transformations (see, for example, [28] ). Alternatively, eqs. (|3.18p and (|3 . 19[) can be 
derived directly form the field equations making use of the Bianchi identities. 



4 The Coupling to Matter Fields 
4.1 The Coupling to the Dirac Spinor 

In this section we study the above theory coupled to a Dirac field. The matter Lagrangian 
of the Dirac field will be taken without any modification with respect to the standard 
one, the only adjustment will be that now, spinorial covariant derivatives will be the most 
general containing torsion, as it was first discussed by Fock and Ivanenko, and as also 
reported for instance in |29j : the Lagrangian is thus given by 



•Cm = 



% — — — 



(4.1) 



where Diip = + ui i Ml/ S^ip and = J^r — ^ ' are the covariant derivatives of 
the Dirac fields, S^ v — g [7^,7^], 7* = 7 M e^ with 7^ denoting Dirac matrices and where 
m is the mass of the Dirac field. From eq. (|4.1j) . it is easily seen that the spin density 
tensor, defined as in section 3, is given by 

V = \$ {7*, Sij} i> = -|fe W Ar (^757^) e£e?e* (4.2) 

In this case, being the spin density tensor (I4.2[) totally antisymmetric, the torsion tensor 
(I3.10[) assumes the explicit expression 

T tsi = nS tsl (4.3) 
where k = l \ k is a constant depending on k (k ^ j). Moreover, from eq. (|3.15p we have 
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By varying (|4.ip with respect to ip, we obtain the Dirac equations 

i-f h D h iP -mip = (4.5) 

Due to eqs. (|4.5[) . it is easily seen that the energy-impulse tensor, derived as in section 3, 
is expressed as 

S y = % - (^D^ - Djftyrf) (4-6) 
The Dirac equations ()4.5[) imply the validity of the identities (see, for example, 30!) 

Vi£« = T^ ik + ^S pqi R p ™ (4.7a) 

W h S ijh = V h S ijh = YP l - E ij (4.7b) 

where the total antisymmetry of the spin and torsion has been systematically used. In 
addition to this, inserting (|4.4[) into (|3 . 1 7[) we obtain the decomposition of the Ricci tensor 
and scalar curvature 

Ra = Rij + - y W ( 4 - 8a ) 

R = R-^-S qpr S qpr (4.8b) 

as well as the identities 

\ (T pqh T^ h gij + 4T ipq T P 3 ")=- {^S pqh S^ h gij - 2n 2 S piq S J qp ^j (4.9a) 

2V p T% = (4.9b) 
2K h fT% = -n 2 S p q S p lq (4.9c) 
In view of eqs. (|4.8[) and (|4.9I) . we can rewrite eqs. (13. 14[) in the form 



\R9ii + ^ 

-kn 2 ^—S pq hS pq>2 gij — 2S P i q Sj' lp ^j + 2nkXt P S P — K^kS p q S p ' iq = S, 

We can decompose eqs. (|4.10j) in their symmetric and antisymmetric parts, that is 

1 - k 2 k 2 /k 2 \ 

~ 2^9ij ^~S pi q S jq p + —S qpr S qpr gij — k I —S pq hS pqh g t j — 2K 2 S piq S J qp \ + 

-K 2 kS p q S P iq = E(y) 



(4.10) 



2 

After some algebraic calculations, eqs. (|4. 1 1[) simplify as 

1 - 



(4.11a) 

K W P S/ + 2nkV p S p tJ = E [<3l (4.11b) 



Rij ~ - 4 W + = E «) ( 4 - 12a ) 

^V P S/ = %] (4.12b) 
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Eqs. (|4.12bl) are clearly identical to eqs. (|4.7b|) , automatically ensured by the Dirac 
equations themselves. This means that the significant part of the Einstein-like equations 
(|4.10[) reduces to the symmetric one (|4.12aj) . The latter can be worked out by giving an 
explicit representation of the covariant derivative of spinor fields. To this end, we denote 
by D the spinorial covariant derivative induced by the Levi-Civita connection. Then 

Ditp = Dii; - ii^ 7 y V (4.13a) 

D4 = D4 + -^,K l]hl h ^ (4.13b) 
Making use of eqs. (|4.2[) . (|4.4p and (|4.13[) . the following identities are easily verified 

ij^Djip - (Djfy jiip = rjrfiDjip - (D^)^ - ^(^7s7 r '0)(V'757rV')5y 

in - ( 4A4a ) 

+ y(V'757i' i /')('0757jV') 

SjS jp h = ~ (^757 r ^) (V>757rV>) 9ij + I (V>757iV>) {jhtHj^) (4.14b) 

S hqp S h ™ = -1 (fa-fif,) (^757r^) (4.14c) 

Inserting the content of eqs. (|4.14p into eqs. (I4.12a[) . we end up with final symmetric 
Einstein-like equations of the form 



1 - i 
2 R9ij = 4 



Tpl(i D ])i> ~ (Dij^hi)^ + (V'757 r '0) (V'TBTt^) 9ij (4.15) 



It is worth noticing that the conservation laws for the energy (|4.7al) reduce to the simpler 

V^W) _ lyi (SnhS* 01 ) = (4.16) 



where E^- := | 



ipjiDjip — (Djip)jiip ■ Comparing with eqs. (|4.14c[) and (|4.15|) . eqs. 
(14. 16)) ensure that the Levi-Civita quadridivergence of the effective energy-impulse tensor 
on the right hand side of (|4.15[) vanishes. 

At the same time, the Dirac equations can be expressed as 

vfbiij} - X'tpj^jgip - imp = (4-17) 

where A := 1 g( 1 3 _ 4fc ) = ff from now on. Eqs. (I4.17[) can straightforwardly be cast after a 
Fierz rearrangement into an equivalent but simpler form given by 

i^ l bii\) + A (^75^75 - V^/>) ip ~ mip = (4-18) 

formally identical to what we would have had if it were in the torsionless case but with 
additional potentials of self-interactions. 

Notice that since in our approach generalizations have place only for the geometrical 
sector, and not for the material sector, then we obtain a generalized form of the gravita- 
tional (metric and torsion) field equations, but not of the matter (Dirac) field equations, 
which remain formally the same we would have had as usual [31] [32]; however, since 
the torsion-spin density field equations are now different, so soon as torsion is replaced 
through these field equations with the spin density of the spinor, the eventual non-linear 
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interaction in the spinorial field equation (the self-interaction of fermions) assumes a form 
that is in general different from that of [21] : although due to the peculiar structure of the 
Dirac spinor these non-linear fermionic interactions in the matter field equations do hap- 
pen to be reduced to the usual form, nevertheless this does not happen for their strength 
as the coupling constant in now different. When in gravity torsion is not neglected, non- 
linearities as those given in [31] appear within the Dirac matter field equations, but when 
the present modification of the torsional action is accounted, those non-linearities get a 
coupling constant that is not necessarily the Newton constant: so it may seem that our 
modification has the only effect of rescaling the torsional coupling constant and nothing 
more. However, in a situation in which one of the most important effects of torsion is its 
induced spin-contact coupling in the matter field equations and in a moment in which tor- 
sion tends to be neglected because such interactions are allegedly small, our modification, 
far from being a mere rescaling of the torsion coupling constant, is what renders those 
interactions relevant at scales in which torsion has always been assumed to be negligible, 
solving the last impeachment torsion would have to face. The possibility to control the 
torsional coupling constant gives one the opportunity to change its sign as well, making the 
torsionally-induced interactions not only important at larger scales, but also attractive. 

In fact, in the present theory as in the Einstein-Cartan theory, the Dirac field equations 
are of the Nambu-Jona-Lasinio type jTD] ; however, in this model the Dirac field has 
self-interactions whose coupling constant can be taken with opposite sign, giving rise to 
Nambu-Jona-Lasinio attractive potentials |33) : in the low-energy limit, this gives the 
possibility to describe superconductivity, therefore providing a geometric interpretation. 

Since our modification affects only the value of the coupling constant, but not the 
structure of the non-linearities, then our model has all advantages but also all disadvan- 
tages of the usual model of Dirac theory with torsion: for instance, the issue of chiral 
anomalies |34) and the fact that the non-linearities are perturbatively non-renormalizablc 
[35) : for these, and maybe others, among all the problems the theory may have we furnish 
no solution more than those discussed in |34[ 135) . or in other works, and therefore our 
modification is neither better, nor worse, than the usual theory discussed in the above- 
mentioned references. By this we do not mean that these problems are not important, 
and on the contrary they ought be analyzed; what we only mean is that these issues are 
untouched by the modification propounded by our work, and so we are not obliged to 
discuss them more than any other paper on torsion for Dirac fields. Our paper is focused 
on the scaling of their coupling constant, and this is the key point we want to deepen next. 

It is important to notice that no matter what this value is, it would at first be thinkable 
to normalize it to the usual value of the Newton constant through the normalization of 
the Dirac field, but a deeper analysis shows that this process would also have the effect 
of changing the scale of the energy density within the gravitational field equations and it 
is therefore unacceptable: so far as our knowledge is concerned, the theory we have here 
constructed is the only one in which the Einstein-Cartan-Dirac theory can be generalized in 
order for the system of field equations, once decomposed in terms of torsionless quantities 
plus torsional contributions written as Dirac field self-interactions, to provide a coupling 
constant for self-interactions of matter leaving gravity unmodified. 

1. Neutrino oscillation in absence of mass states. The first example we will 
consider in the case of a coupled system of fermions is given by the simplest couple of 
spinors, that is the pair of semispinors both massless; these semispinors without mass will 
possess spin-torsion coupling providing some sort of mixing between them despite the fact 
that in the accepted model neutrinos only have massive oscillations: the aim is therefore 
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to compare this case with the accepted model to appreciate their discrepancies. 

Considering the field equation for two scmispinors without mass v\ and v 2 we have 
that their field equations are given by the usual field equations for semispinors taken to 
be left-handed and so in the purely massless configuration 



(4.19) 
(4.20) 



in terms of the most general covariant derivatives that can be decomposed into the 
simplest covariant derivatives as in the torsionless case plus torsional contributions as 



iYD^vi + Xv^pV^vx = 
ij^D^ + \v x 7^17^2 = 



(4.21) 
(4.22) 



with spinorial interactions of each semispinor with the other; these have now to be written 
in an alternative form by employing Fierz identities V'7 M V'7 A1 X = —i'l^Xl^i' valid for any 
couple of left-handed semispinors ip and \ identically: after this Fierz rearrangement we 
get the previous field equations transcribed as 

i^D^i - | (z/i7^i - ^27^2) 7^1 - x (^27^1) 7^2 = (4.23) 
i^D^vz - ^ (^17^2) 7a^i + f ~ ^27^2) 7^2 = (4.24) 

in which the interactions of each semispinor with the other has been written in a form 
that is particularly interesting; in fact, such a coupled system of field equations is formally 



^7' J 



D, 



A ( (fl7fi"l ~ ^27^2) 
2(2^17^2) 



2(^17^2)* 



= 



(4.25) 



in which we achieve a considerable compactification. Next we simply have to call the 
doublet of semispinors 



with which to define the triplet of vectors 



\v^^av = A p 



so that we have 



or equivalently 



L ^ + iX {A] A }iAl 



A/j, iA^ 



=0 



D^v + iXA^ ■ av 



= 



as it can be checked; then upon introduction of 

= b^v + 2i\A^ 

we finally have 

i-y^v = 



(4.26) 



(4.27) 



(4.28) 



(4.29) 



(4.30) 



(4.31) 
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in which the free covariant derivative is now part of a gauge covariant derivative U^, 
and thus this is the field equation in the torsionless case but in presence of an additional 
gauge interaction. To see that this is actually a gauge invariant interaction, notice that as 
the doublet of semispinors v transforms according to the St/ (2) group 



(4.32) 



then the triplet of vectors and the triplet of Pauli matrices that generate the infinites- 
imal transformation of the SU(2) group are combined into ■ a transforming as 



A„ 



A a - d u 



-2iX8- 



(4.33) 



as the gauge connection of the SU(2) group and field equations (|4.3ip are invariant for 
SU (2) local transformations [36]. Here A is to be tuned to the neutrinos' oscillation length. 

The field equations (|4.3ip are still the field equations for semispinors in the massless 
configuration but now due to the masslessness of the neutrinos their covariant derivative 
contains interactions gauging the SU(2) group: so that the two types of neutrinos will 
be converted into one another according to a sort of massless oscillation [36 . This is 
a remarkable discrepancy between the present theory, in which it is precisely because 
neutrinos are massless that they may oscillate, and the commonly accepted theory, in 
which only massive neutrinos can possibly oscillate. 



2. Leptons with weak interactions and Higgs field. Our next example will be 
the immediately more complex one, that is the coupled system of fermions, one of which 
is a spinor with mass while the other is a semispinor thus massless; these spinors with 
massive states will possess spin-torsion coupling resulting in interactions between the two 
fundamental fields while from the point of view of the standard model the couple of electron 
and neutrino only has the weak interactions: the aim is therefore to compare this case 
with that of the weak forces and see whether or not similarities may arise. 

In this case the field equations for the electron and neutrino fields e and v are given 
by the previous field equations, but because now the system is constituted by a couple of 
spinors then the spin will be the sum of the two spins and the spinorial field equations are 

i^^D^e — A (e7 M e7 M e + Vj^i/^^e) — me = (4.34) 
vfDpV - Ae 7p 7 5 e7'V = (4.35) 

in which the fact that the electron is massive while the neutrino is massless is the reason 
that prevents these two fields to mix into a doublet; after by employing a Fierz rearrange- 
ment as we have done before, they can be cast into the form 



+q tan eZ^e - j^Z^e L + -%W*^v -He-me = (4.36) 



vy^Dpe + 2A(cos#) 2 e7 5 e7 5 e 

n"D^ +J ^Z^ + ^W,reL=0 (4.37) 

once we define 

Z» = -A [2(sin(9) 2 e7^e - e Ll »e L + v^v\ (4.38) 

4(sinfl) 
q\f2 sin 



W* = -\{e Ll » V ) (4.39) 
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and 



H = Aee2(cos6») 2 (4.40) 

showing that field equations (|4.36p and (|4.37[) are formally identical to the system of 
field equations for the lepton fields after the symmetry breaking in the standard model, 
although both weak and Higgs boson fields are here composite [37]. Notice that because 
the field equations of the leptonic fields are known, it is possible to compute the divergence 
of the leptonic currents (|4. 38114. 39"]) . for which at the weak interaction scales a reasonable 
estimate gives the partially conserved axial currents 

[(i+^)Vu^+2^v„*q = -2s*a? (4.4i) 

71 K 1 + £ )(V^^+i g tan^Wi,)+2WV„f ] = ^ (4.42) 

which in the limit of the approximation are precisely the partially conserved axial currents 
one would have had in the standard model so soon as the torsional coupling constant is 
tuned as 4A = 4j in terms of the Higgs vacuum [35]. Notice that with this fine-tuning 
for the constant, the weak interactions are reproduced not only in structure but also in 
strength while the masses of the composite mediators have not only the proper ratio but 
also the correct value. This would imply that among theories coming from torsion and the 
standard model, so long as we take into account leptonic weak scattering and the weak 
mediators masses it is not possible to appreciate any discrepancies, and we have to go 
to higher energies to see that in our model the weak mediators must display an internal 
structure that can never be present for the standard model weak bosons [39 . 

Similar results have also been obtained both for leptons and hadrons in circumstances 
in which the underlying symmetry was still unbroken [401 141) . 

3. Dirac field as a Van der Waals gas. We will now leave the treatment of the 
coupled system of Dirac fields to focus on the single Dirac field; writing the Dirac field 
equation in the standard representation allows us to obtain the slow-speed weak-field 
approximation 

*ff + ^<r k Dk<T a D a cj> - A (0t0) <p - m4> = (4.43) 

representing cold-matter fields with self-interactions with the structure of an energy term. 

Thus in the present theory we have that the Dirac field equations have the form of the 
Nambu-Jona-Lasinio field equation further approximated to the Pauli-Schrddinger field 
equation of the Ginzburg-Landau type [TU]. In it the coupling constant A is to be set on 
the value of the specific condensed state system we would eventually like to study [3"3"] . 

These field equations have energy levels that can easily be computed to be 

E = m + ^ + Xu 2 (4.44) 

where u 2 is the square of the field, representing the density of the field itself, or equivalently 
in the quantum-mechanical interpretation the probability to find a give particle or again 
for collective states the number of particles per unit volume u 2 = on the other hand, 
the usual thermodynamic interpretation of the kinetic energy is the temperature T of the 
gas under consideration: considering both these interpretations, the energy levels (I4.44p 
receive the quantum-mechanical thermodynamic interpretation of 

E = T+X^ (4.45) 
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which is the well known expression for the energy level of the Van der Waals gas. 

By employing the equation of the energy (§^) T = (jpp) v T — P it is easy to see that 
the equation of state gives the pressure P as a function of volume and temperature as 

(P-\$)F(V) = T (4.46) 

as the Van der Waals equation of state with generalized volume factor F(V) and with 
coupling constant A determining the type of interaction of the gas: for the commonly 
accepted positive values the Van der Waals pressure is positive resulting into a repulsion, 
while for negative values the Van der Waals pressure is negative resulting into an attraction. 

The interpretation of a Dirac field or a collective state of Dirac particles as a Van der 
Waals gas can be read by interpreting the torsionally-induced self-interactions as the Van 
der Waals interactions among different particles of the gas, as it might have been expected 
from the fact that both type of interactions drop with a \ dependence . 



4.2 The Coupling to the Fluid with Spin 

The importance of the Dirac field is mainly due to the fact that it is the simplest quantum 
field that can be used in the description of natural phenomena, including even cosmological 
applications (see for instance [13]); other spin-i spinor fields like ELKO may also be used 
in cosmology such as those tackling the problem of Dark Matter [16] . In the present paper 
however, we wish to employ yet another matter field to describe large-scale structures: the 
fluid with spin, that is the Weyssenhoff fluid [H] |3S] . 

The spin fluid is defined by its energy-momentum tensor, of the form 

= WP J + p (WW - g l] ) (4.47a) 

and its spin density tensor given by 

S tJ h = S lJ U h (4.47b) 

where W (U l W = 1) and P l denote respectively the 4-velocity and the 4-vector density 
of energy-momentum, while S%j — —Sji is the spin density of the fluid: the 4-velocity and 
the spin satisfy the convective condition 

SijW = 0. (4.48) 

Making use of eqs. (|3.15l) . (|4.47b|) and (|4.48l) . we obtain the expression of the contorsion 

K ui = i [-A(k) (S sti + S Us + S ist ) + B(k) {-S tsi + S slt - So,)] (4.49) 

where 

— 2k 1 

= 7 ttt rr and B(k) = -, (4.50) 

v ' (l-4fc)(l + 2fc) v ; (l + 2fc) y ' 

and due to the convective condition (I4.48[) . it is seen that torsion and contorsion tensors 
are both traceless. From this, we can write the Einstein-like equations in the form 

Rii ~ \R9ii - Vh*/ - KjK Jp h + l -K h ^K lp h g l3 + 

k ( 4 - 51 ) 
- 2 {T mh T^ h 9ij + m pq T P 3 q ) + 2kV h T% + 2kK h fT% = 
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It is easy to verify that the antisymmetric part of (|4.51[) 

- V h K [Tl] h + 2kV h T h m = E M (4-52) 

amounts to the conservation laws for the spin (13.191) . Indeed, in view of eqs. (|3.16|) and 
(|4.49[) . eqs. (|3.19[) are seen to reduce to 

iv^ + Efc] =0 (4.53) 

Inserting eqs. (|3.5I) into (|4.53[) and taking the identity = — K^ h into account, 
we end up with eqs. (14.521) . 

Saturating eqs. (|4.53[) with Ui we obtain the explicit expression for the vector density 
of energy-momentum 

P 3 = pU j + S ij U h V h Ui (4.54) 

where p :— U % Pi. Inserting eqs. (|4.54j) into (|4.47ap . we get the form of the energy- 
momentum tensor 

Y, ij =(p+p) WW - pg ij + WS pj U h V h U p (4.55) 
The significant part of the Einstein-like equations turns out to be the symmetric one 

Rii - \R3ii = E W ) - ^l\ k) li ) +2 k) SpqSPqUiUj+ , , 

(l-6k) 1 (4 ' 56) 

;S Pq S^ 9lJ - — V fe (S\U 3 + S^Ui) 



8(l-4fc)(l + 2fc) pq * lJ 2(1 + 2k) 



4.2.1 FLRW cosmological models 



In the above described theory, Friedmann-Lemaitre- Robertson- Walker (FLRW) cosmolog- 
ical models 

ds 2 = dt 2 ( dr 2 + r 2 M 2 + r 2 sin 2 q ^2) ^ ^ 



(with K = —1,0,1) can be considered under the hypothesis that the isotropic and homo- 
geneous universe is filled with an unpolarized spinning cosmological fluid. Being the spin 
randomly oriented, we can assume that the average of the spin and its gradient vanish, 
but the same is not true for the spin-squared terms < S pq S pq >.— \S 2 . The conclusion 
follows that, after averaging, eqs. (I4.56[) reduce to 



(!-6*0 g SPq 



(4.58) 



8(l-4fc)(l + 2Jfe) 
and the energy-momentum tensor is 

£« = {p + P) UiUj - p gij (4.59) 
The Friedmann equations derived from (|4.57|) . (|4.58p and (|4.59|) are 

a = - \-p - 3p + 4C(k)S 2 ] (4.60a) 
6 
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a 

T 



^ + P -C(k)S 2 



(4.60b) 



where we have defined C(k) :— 4 ^ J^^i^fc) • Under the stated conditions, it is easily seen 
that the conservation laws (|3.18p and (|3.19l) give rise respectively to the equations 



p + 3(p + p) 



(4.61a) 



S + 3S- 







(4.61b) 



Supposing again a state equation of the kind p = \p (0 < A < 1), eqs. (|4.6ip admit general 
solutions 



p a 



-3(1+A) 



S 



(4.62) 



Po and So being suitable integration constants. Eqs. (|4.61l) ensure that the quadridiver- 
gence (with respect to the Levi-Civita covariant derivative) of the Einstein-like equations 
(|4.58|) vanishes, while eqs. (|4.62[) yield the relation 



S 2 = S A 



(4.63) 



We can see from the previous equations that the spin contributions introduce two impor- 
tant modifications in cosmology. The first one is that torsion may open closed universes, 
as a matter of fact from equation (|4.60b[) it follows that the critical density is defined by 



p c = H 2 + C(k)S 2 
accordingly, those cosmological models with density p such that 

H 2 <p<H 2 +C(k)S 2 



(4.64) 



(4.65) 



are closed in general relativity, but open in our theory. 

The second important modification is that from eq. (|4.60bl) is easily seen that the spin 
contributions can avoid the initial singularity. Indeed, supposing for simplicity the metric 
(|437|) spatially flat (K = 0) as well as C(k) > 0, from eqs. (|4T60b|) . gib]) and fi^g]) we 
derive the existence of a minimum value for the scale factor 



C{k)S 2 ^ 3 ^- y > 



Po 



corresponding to a maximum value pi of the energy density 



(4.66) 



Po 



C(k)S 2 



1 + A 
1 — A 



(4.67) 



Moreover, from eq. (|4.60a[) it is seen that we can have an initial accelerated expansion 
(a > 0) for values of energy density 



p> p f 



" (1+3A) ^ 
4C(k)S 2Po 



l + A 
1-A 



(4.68) 
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The value of the scale factor a/ associated with pf is 



Since pi > pf, the accelerated expansion stops when density pf and factor scale a/ are 
reached. It is worth noticing that, while the values of aj and a/ (as well as of pi and pf) 
depend on the parameter k, their ratio afja% does not and it is identical to the result 
already obtained in the standard Einstein-Cartan framework |46) . 

In view of this, the cases K ^ seem then more interesting. To discuss this point, we 
limit to consider a fluid with an equation of state of radiation kind A = 1/3. Accordingly, 
from eq. (|4.60b|) we find the minimal value Oj by solving the equation (d = 0) 

C(k)S^-p ^+3K = (4.70) 



The solutions are 



1 _ Po ±y/p>-12KC(k)S* I 2C(k)Sj 

2C(fc)5 2 \j p ±^/pl-12KG{k)Sl { ' ' 

When K = 1, these solutions make sense only if 

C(k) < ^ (4-72) 



125, 



representing a further restriction on the parameter k (together with the requirement 
C(k) > 0). The two solutions (|4.7ip coincide with the minimum and the maximum 
reachable values of the scale factor of the universe. 
For K = — 1 the imaginary solution 



V po - y/p 2 + \2C(k)Sl 
has no physical meaning. So there is only a minimum 



2C ^ (4.73) 



2C(fc)S'n 

(4.74) 



po + ^pl + \2C(k)Sl 



Anyway, in both cases K = ±1, from eq. (|4.60a[) it is seen that the universe can undergo 
an accelerated expansion for a; < a < a/ , where 



In the case of a closed universe the final expansion is limited by the condition (|4.72[) . In 
the case of an open universe the ratio 



a/ = / po + y/^ + 12C(fc)^ 
a 4 V Po 
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depends actually on the parameter k. Note that C{k) — > +00 for k — > — \ + . We conclude 
that in this case the parameter k may be chosen in such a way that the accelerated 
expansion is sufficiently long to predict a flat observed universe. 

At the end of this section with add some final remarks about the horizon problem, one 
of the classical problems of cosmology. 

It is well-known that in the general relativistic cosmological models the horizon problem 
is a consequence of the initial singularity. In fact, in general relativity, from the Friedmann 
equation 

3K 



a 2 



(4.77) 



it follows that a — > 00 when a — » 0. This velocity must be compared with 



d i = ±c ~ ^ 

at a 

which has been obtained by imposing the condition ds 2 = 0; this means that the signals 
exchanged in a fluid with equation of state p = Xp by different parts of the universe cannot 
go faster than c/a , while the various regions of the universe move away more rapidly, when 
a is very small. It follows that the particle horizon 

l(to)=a(t ) f° ^- (4.79) 
Jo a{t) 

is finite, while H ~ 1/t is not an integrable function in the same interval. The consequence 
is that there are regions of the universe which are not connected causally. In General 
Relativity this problem has been solved, together with the flatness problem and other 
classical problems by introducing the inflationary scenario. 

But comparing equations (I4.77|) and (|4.78|) we see that as the universe expands, there 
is an epoch when 

->v = -d (4.80) 

a a 

where v is the recession velocity related to the distance d by the Hubble law. So assuming 
in our theory a similar expansion law, but without the initial singularity and assuming 
that C(k) is large enough, we can show that horizon defined by equation (|4.79[) can cover 
a region large enough such that all the observed parts of the universe are in mutual causal 
connection. This explains why all the regions of the universe show the same temperature 
and the same physical properties. 



5 Conclusion 

In this paper, we have constructed a generalization of the Einstein-Cartan gravity in which 
both curvature and torsion are considered each with its own coupling constant, the one 
related to the curvature being the gravitational constant while the one related to torsion 
still undetermined; we have applied this geometrical background to the case of Dirac fields 
and spin fluids, showing that the corresponding equations have torsional contributions 
whose coupling constant can be suitably chosen: in the case of the coupling to Dirac 
spinors, we have shown, through two examples, that the torsional contributions with a 
properly tuned coupling constant make the non-linearities of the Dirac equation relevant 
already at subatomic scales; in the case of the spin fluid, where in cosmological models the 
torsional contributions were already present, the tuning of the coupling constant amplifies 
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them even more: in particular, in the case of spatially hyperbolic FLRW model, the 
accelerated expansion lasts longer, solving the horizon problem in a natural way. 

A problem this approach may not be able to address is the fact that there does not 
seem to be a unique value of the constant for which all of the applications we mentioned 
fit into the presented frame, and so only an even more general theory of gravity may be 
able to provide a running coupling constant, possibly scaling with the energy, such that 
all applications above, and maybe more, can fit into a single scheme |30| . 

If this approach of ours, as well as any of its generalizations, really works, then torsion 
would not only be observable, but it might have already been observed, although we have 
not been able to recognize it for what it is. 

We are aware of the fact that torsion may not be the answer to all, and maybe not 
even to most, of the problems of physics, but it may possibly be for some of them, and 
our theory can allow us to see how. 
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